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Electronic states at domain walls in bilayer graphene are studied by analyzing their four and two 
band continuum models, by performing numerical calculations on the lattice, and by using quantum 
geometric arguments. The continuum theories explain the distinct electronic properties of boundary 
modes localized near domain walls formed by interlayer electric field reversal, by interlayer stacking 
reversal, and by simultaneous reversal of both quantities. Boundary mode properties are related to 
topological transitions and gap closures which occur in the bulk Hamiltonian parameter space. The 
important role played by intervalley coupling effects not directly captured by the continuum model 
is addressed using lattice calculations for specific domain wall structures. 
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The electronic properties of few layer graphene systems 
depend sensitively on the atomic registry between neigh- 
boring layers [T]. For bilayer graphene (BLG) with AB 
stacking interlayer hybridization of orbitals on eclipsed 
lattice sites gaps out half of the low energy degrees of 
freedom, replacing the pseudorelativistic description of 
single layer graphene by a low energy theory in which 
quadratically dispersing chiral bands touch at discrete 
points in momentum space ^2 . A perpendicular elec- 
tric field further breaks inversion symmetry and creates 
a semiconductor in which the gap size is determined by 
the electric field magnitude [2]-[5] and saturated at the 
strength of interlayer hybridization The possibility 
of exploiting this type of field tunable gap is being vig- 
orously pursued in ultra-clean dual gated devices [6Hl0]. 

More recently it has been appreciated that the field- 
induced gap admits a topological interpretation I12j . 
The low energy theory for BLG can be represented 
by an effective two-band model from which it is read- 
ily seen that inversion-symmetry breaking induces large 
momentum-space Berry curvatures [121 113j . The Berry 
curvatures have opposite sign near the two inequivalent 
Brillouin-zone corners (valleys) at which gap is opened, so 
the integral of the Berry curvature over the full Brillouin 
zone is zero. Nonetheless, the integral of the Berry cur- 
vature within a single valley is nonzero and this allows 
a topological analysis of the valley-projected electronic 
spectrum. This idea has been developed in a continuum 
analysis of the subgap electronic states bound to a BLG 
domain formed by a sign reversal of the electric field be- 
tween the layers [TTJ [T4lll7j . These electric-field walls 
(EFW) are predicted to bind pairs of subgap chiral co- 
propagating boundary modes, an interesting feature that 
can be related to the change in sign across a domain wall 
of a valley-projected topological index. 

In this paper we examine a related BLG domain wall 
problem in which the interlayer electric field is uniform 
but the layer stacking switches from AB to BA. This new 
version of the problem changes the boundary conditions 
for matching the electronic states of the two bounding 



phases and requires that we augment the two-band model 
of BLG (51[TTJ[T2] by accounting for all four sublattice de- 
grees of freedom. Nonetheless we find that layer-stacking 
walls (LSW) bind electronic states with the same chi- 
ral structure as for the EFW studied previously. We 
make this connection explicit by mapping the two prob- 
lems onto each other within a family of four-band BLG 
Hamiltonians. Our results demonstrate that the topo- 
logical transition in a LSW structure is associated with 
a finite momentum gap closure in the parameter space 
of four-band BLG Hamiltonians. We construct a phase 
diagram (Fig. [I]) which identifies the different types of 
topologically protected states that appear in BLG sam- 
ples in which both the interlayer electric field and the 
layer stacking order vary in space. This analysis iden- 
tifies yet a third type of domain wall in which the two 
pairs of chiral modes within a single valley are coupled, 
gapping the spectrum and annihilating boundary modes. 
Our results are supported by a continuum analysis of the 
domain wall states, lattice calculations for specific de- 
fect structures, and analysis using quantum geometrical 
arguments. Taken together these elements provide a gen- 
eral framework for understanding the origin of the valley- 
projected topological states in BLG, and their fragility in 
the presence of intervalley scattering. 

The electronic states for BLG can be represented 
by four-component wavefunctions ^ = {iPatj4'Btt''Pabi 
4'Bb ) where -0 denotes the atomic orbital centered on the 
A B sites of the top or bottom layer. At low energies 
the Hamiltonian can be expanded for small q around the 
two inequivalent Brillouin zone corners: HivK + q) with 
u — ±1 denoting K and K' . Using Pauli matrices cr 
to represent operators that act on the sublattice degree 
of freedom within a layer and x to represent operators 
acting on the layer degree of freedom, the BLG Hamil- 
tonian can separated into layer diagonal and layer off- 
diagonal contributions by writing T-L = Hq + Hint- We 
find that for AB stacked BLG in which the A sites of 
top layer hybridize with the B sites of the bottom layer. 
Ho = vq^a^TQ + qyUyTf) and "Hint = -l{crxTx - cfyTy)l2 
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FIG. 1. A phase diagram illustrating the distinct valley- 
projected topological phases of BLG and the critical lines that 
separate them. A sign change of A denotes a reversal of the 
interlayer electric field. A sign change of fi denotes a tran- 
sition from AB to BA interlayer registry. The spectrum is 
gapped except along the lines A = and /j, = 0. 



where energies are in units oi hv, 7 = 'ji/hv, 71 is the 
near-neighbor interlayer hopping amplitude, and v is the 
electron velocity in an isolated layer. For the reversed 
BA stacking order the interlayer coupling term becomes 
Hint — l{<^xTx + f^y'''y)l'^- When present an electric po- 
tential difference V adds ActoT2 with A = V/2hvp to the 
Hamiltonian. 

When 7 ^ A it is convenient to eliminate the high 
energy degrees of freedom at ±7 to arrive at an effective 
low energy two-band model [5] 



'Hy =gy{q)-cr, 



(1) 



where the a matrices act on two component spinors 
{iPbtt^Ab) ill tli6 low energy subspace for AB BLG and 
a^ig) = {-{qI - <ll)h,'^mxqyh^^)- Eq. ^ admits a 
geometrical interpretation in which the negative energy 
eigenstates are spinors aligned with —Quio) and the filled 
band has a momentum space Berry curvature \12\ I13j 
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Because of the v dependence in Eq. ^ the integral of 
57i/(q) over the full Brillouin zone is zero and the filled 
valence band carries total Chern number = as re- 
quired by time reversal symmetry. However, for small A 
the Berry curvature is strongly peaked at the gap min- 
ima near K and K' . Consequently, the integral of ^^{q) 
over an individual valley is accurately defined and the 
valley Chern number = — j/sgn(A) = ±1. The valley 
Chern number changes by AN^, = ±2 across an EFW 
which can be associated with the appearance of pairs 



of valley-projected edge modes co-propagating along the 
boundary. These chiral modes have been obtained by 
analytic solution of the low energy two-band model in 
the presence of a sharp EFW and by numerical solution 
for a spatially varying A(r) that smoothly connects two 
electric-field reversed states [TT]. As noted in previous 
work [12 E] , the introduction of a valley Chern number 
in this context is approximate since strictly speaking the 
construction does not map the full periodic Brillouin zone 
onto the parameter space of T-L^. Nonetheless, when A is 
small and intervalley scattering is absent the computed 
change AN^ can be interpreted as a topological quantity, 
since fl,^{q) is integrated over a closed surface produced 
by "gluing together" two integrals for the individual N^, 
along a common boundary. 

We now turn to the case of a LSW at which the bi- 
layer registry reverses from local AB to local BA with 
A held constant. Crossing a LSW changes the interlayer 
coupling matrix Hint and switches the orbitals that span 
its low energy subspace. In this case evaluation of the 
Berry curvature requires consideration of all four degrees 
of freedom in the bilayer Dirac problem. Alternatively, 
one can identity the topological origin of LSW modes 
by examining the residual phase twists induced at large 
momentum q in the eigenstates of the generalized Hamil- 
tonian, 
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Hlsw = At^ + vq^a^ + qyOy + - {a^r^ - fiiTyTy) ,(3) 

which reduces to the AB (BA) forms when fi 1 (— 1). 
For q I A 1, 7 degenerate single layer states '^^„{q) = 
{'4'ny,T, i^fj.iy,B) deep in the filled band with energies E = 
— \q\ are split by A and are mixed by 7 in the projected 
Hamiltonian 



|g|Ao + AA^ 



(e*'''"*A+ + e-'^'^'^A,) ,(4) 



where A are 2x2 Pauli matrices acting in the 4"^!, sub- 
space and (j) = a,Tctaii{qy / qx) ■ For large q the eigenstates 
^ij,v,± written in the original four orbital basis are 
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where ^ = ^7^/4 -I- A^ and we explicitly display the 
overall U{1) phases a^y,±. Using Eq. ([5| we calculate 
the momentum space Berry connection 
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The change in the valley Chern numbers upon passing 
from the /i = lto/i = — 1 states is obtained from the 
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FIG. 2. Propagating chiral boundary modes (red) near the 
K (top) and K' (bottom) points calculated using the four- 
band continuum model for a domain wall separating regions 
with local AB and BA stacking registry in the presence of 
a uniform layer-potential difference A. The shaded regions 
show the support of the continuous bulk spectrum as a func- 
tion of momentum q parallel to the domain wall. The results 
were obtained for the parameter values A = 0.05, 7 = 0.1. 



FIG. 3. Top: Comparison of valley K domain wall spectra 
for electric-field and layer-stacking walls. Both walls support 
a pair of co-propagating chiral modes. Bottom: Interface 
spectrum calculated with a four-band theory for a topologi- 
cally compensated boundary at which both layer stacking and 
electric field change sign. In this case the boundary spectrum 
is completely gapped even in a continuum model. 



loop integral of the trace of 
indices ±, which reads 



-/i^^i/^± over jjL and band 



2v + m^i^y^+ + m_i_,,^_ - mi^y^+ - mi^^^_ (7) 



Here m^i,,± are integer valued winding numbers of the 
overall phases a^^,±. The A dependence of this re- 
sult vanishes after tracing over the filled bands, demon- 
strating that the valley Chern number in BLG is shared 
among all the occupied bands rather than being confined 
just to its low energy states as is often assumed. AiV^ is 
a topological index provided that the difference is eval- 
uated in the same gauge for the two bounding phases, 
which requires that m^i,^^± — mi^y^± for this boundary. 
It follows that AA^^ = 2v and therefore that a domain 
wall separating insulating regions with local AB and BA 
registry will also confine pairs of valley-projected chiral 
modes propagating along the boundary with opposite ve- 
locities in the two valleys. Fig. [5] confirms this result 
by showing the spectra calculated by matching the full 
four component wavefunctions of Eq. ([s]) across a sharp 
boundary where /i switches from 1 to —1. 

Using Eq. ([5]) we find that at large q the wavefunc- 
tions on the two sides of the LSW are related by a gauge 
transformation, 



(8) 



with a different phase twist induced in each layer. It 
follows from the accumulation of internal phases in the 
two layers that AiY^, = 2^v. EFW walls, where the sign 
of the potential difference between layers A switches but 
the atomic registry does not change, can be analyzed 



similarly. Although A does not appear explicitly in AiVj^ 
in Eq. ([t]) there is an implicit dependence through the 
C/(l) phase prefactors. We find 



■Hl5w(-A,/x) = T.J;'HLSw{^,-^J')Tx 



(9) 



i.e. that a sign reversal in A can be absorbed in a sign 
change of /i combined with a change of basis. Using this 
construction the negative energy eigenstates at large q 
on the two sides of the EFW are related by the following 
[/(I) gauge transformation 

^if^^^ ^ ^^e'M-'^^^^i , (10) 

which produces the same overall AN^, = 2^,v for the 
EFW. Eqs. ([8|-([l0| compactly express the relation be- 



tween these two different types of domain wall in the 
four-band theory. This is also illustrated in Fig. |3] which 
compares the valley K spectra computed for an EFW and 
a LSW showing their common chiral boundary modes. 

In the LSW case, unlike the EFW case, analyzing the 
continuity of wavefunctions across the interface requires 
consideration of all four bands. The common topologi- 
cal origin of the domain wall spectra therefore becomes 
apparent only in a four-band continuum theory. Nev- 
ertheless, by integrating out the high energy bands at 
energies ^ ±7 we are able to construct a two-band ef- 
fective model away from the domain wall in which for 
either case AA^^ is assigned to a sign change of the Berry 
curvature of the lower band. In this approach, Eq. ([T]) 
reads g^{q) = {-{ql - q^)/!, ^^ivq^qy/^, A) with or layer 
Pauli matrices that act on different spinors in the /i = ±1 
cases: on {iI^Bt:4'Ab) for = 1 and on (-0^^ , ) for 
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/i = — 1. Because of inversion symmetry breaking, the va- 
lence band acquires a momentum space Berry curvature 

mm, 



(q4+ ^2^2)3/2 ' 



(11) 



which integrates over a single valley to — /ii'Sgn(A). Ob- 
viously, the valley Chern number changes by two across 
either a EFW or a LSW. Based on the bulk-boundary cor- 
respondence, pairs of valley-projected edge modes should 
co-propagate along the interface. 

The A — plane phase diagram in Fig. [T] identi- 
fies distinct BLG topological phases. Phase bound- 
aries occur along the ^ and A axes where the spec- 
trum of 'Hlsw{^t undergoes gap closures at = 0. 
'Hlsw{^ = 0,^ = 1) describes an ungapped BLG sys- 
tem in which quadratic band crossing occurs exactly at 
q = 0, as seen in the left panel of Fig. |4] The gaps that 
open for the case of A 7^ and /i = ±1 are the elec- 
tric field induced gaps easily understood within a two- 
band model. The boundary with A ^ and ^jl — Q also 
has a gap closure, but it occurs at two finite momenta 
qx = ±-\/A2 -I- 7^/4 along the Qy = line where band 
crossing is possible because ax is a constant of the mo- 
tion. For deviations in either or qy the degeneracies at 
the band-crossing points are lifted at linear order, imply- 
ing the conical gap closure illustrated in the right panel of 
Fig. |4j When /i 7^ ±1 the original quadratic gap closure 
fissions into a pair of linear Dirac singularities each of 
which carries half the original winding number. Trajec- 
tories in the Hamiltonian parameter space that connect 
these topologically distinct ground states and involve dif- 
ferent parameter values can shift the momenta at which 
the gap closures occur, but cannot eliminate them. For 
example, when A 7^ but 7 = in Eq. ([3| the layers 
decouple and the gap closure at £■ = degenerates to a 
closed Fermi ring with radius A. 

Fig. [1] also illustrates the possibility of a third type 
of compensated domain structure (labeled c) at which 
both the layer registry and interlayer electric field are 
reversed. Variation of local band parameters along this 
line connects two bilayer states that are distinct but have 
AN^, = 0. As illustrated in the the lower panel of Fig. |3j 
spectra obtained by matching solutions across this com- 
pensated domain wall demonstrate that it hosts two pairs 
counter-propagating modes within the same valley, which 
hybridize and completely gap the spectrum. 

The continuum model is able to explain the topolog- 
ical origin of the gapless interface modes. However, the 
short-range physics near the domain wall, which may be 
of essential significance, is not captured in the contin- 
uum Hamiltonian. Importantly, the single valley physics 
that protects the chiral domain wall solutions can be pre- 
empted by sufficiently strong large momentum scattering 
that acts to recouple states in the two valleys. In fact, 
Fig.[2]suggests that these single valley domain wall modes 
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FIG. 4. Bulk band structures of three forms of BLG. Left: 
BLG with uniform AB registry and no interlayer electric field 
characterized by quadratic touching between its two low en- 
ergy bands. Middle: the field-induced gapped phase in which 
the BLG degeneracy is lifted by an interlayer electric field. 
Right: a critical state with nonzero interlayer electric field 
where gap closure occurs at ^ = as BLG crosses from AB 
to BA registry. The zero energy gap closures in the left and 
right panels separate four distinct gapped phases of BLG as 
shown in Fig. [l] 



ultimately reconnect with each other. To study this fur- 
ther we construct a specific lattice model and use it to 
investigate how both lattice and interfacial effects, which 
couple the two valleys, influence the domain wall modes. 

As depicted in Fig. [5] we consider the simplest LSW, 
i.e., a grain boundary separating BLG into left and right 
domains. Near the LSW, the lattices are continuous in 
one layer but fractured along a zigzag edge in the other. 
This introduces additional zigzag boundaries in the bro- 
ken layer and allows switching of the bulk stacking order 
from BtAb (/i = —1) on the left to AtBb (/i — 1) on 
the right. For comparison, we first calculate the band 
structures for the case of uniform gapped BLG and for 
the case of gapped BLG with an EFW at which stacking 
order is preserved. As expected and shown in Fig. [6|^a) 
and (b), quantum valley Hall edge states [T^] and two 




FIG. 5. The simplest LSW separating BLG into a left (L) 
domain with BtAb stacking (/^ = — 1) and a right (R) domain 
with AtBb stacking (/i = 1). When the BLG is uniformly 
gapped, gapless modes emerge along the outer zigzag edges 
(E) as well as along the LSW interfaces (I). The lattices are 
continuous in the bottom (B) layer but have a straight crack 
in the top (T) layer. The dashed lines denote the tunneling 
between the domains within the continuous layer. 
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FIG. 6. Gapless modes for (a) uniform gapped BLG, for (b) 
gapped BLG with an EFW, and for (c) gapped BLG with a 
LSW as depicted in Fig. [5] The yellow states localize on the 
outer zigzag boundary and they are doubly degenerate in (b) 
and (c). The green states localize on the EFW in (b) and 
on the LSW in (c). In (c) the green (magenta) LSW states 
localize on the broken (continuous) layer. To illustrate these 
different cases, we choose parameter values: t = 1, 71 = 0.3, 
and F/2 = 0.25. 



flat bands appear at the outer zigzag edges in uniformly 
gapped BLG. In the sample with an EFW there is an 
additional pair of co-propagating chiral gapless modes 
which emerge at each vaUey. Fig. ^c) shows the situa- 
tion for a LSW with a uniform interlayer electric field; 
surprisingly there are three instead of two gapless modes 
per valley in this case. 

We investigate this problem further by studying the de- 
pendence on the tunneling amplitude tc across the LSW 
shown as the dashed lines in Fig. [5] Without tunnel- 
ing (Fig. [7][a)), the boundary mode spectrum yields two 
copies of the gapped BLG spectrum shown in Fig. |6][a), 
and thus there are two chiral gapless modes in each valley 
as anticipated by the continuum model. The flat bands 
represent the states localized on the grain boundary lines 
B^^ and A^^ in Fig. [5] The leading effect of turning on 
the tunneling is that the pair of degenerate flat bands 
(magenta bands in Fig. [7| are split and become disper- 




FIG. 7. Evolution of gapless modes for gapped BLG as 
a function of hybridization across the LSW. The tunneling 
parameter tc is respectively in (a), 0.2 in (b), 0.6 in (c), and 
1 in Fig. |6|c). The other parameters have the same values 
as in Fig.Tm In panel (a) all the yellow, green, and magenta 
states are doubly degenerate. 



sive, as described in Fig. and (c) . When the tunnel- 
ing is larger than the electric field induced gap, the flat 
band split downward is pushed down to the valence band 
and becomes the third gapless mode shown in Fig. [6]Jc). 

The other two gapless modes (green bands in Fig. [?]) 
localized on the grain boundary of the broken layer are al- 
most degenerate due to the inversion symmetry between 
the lines Bj^^ and Aj;^ in Fig. [5] This degeneracy is ex- 
act at A;a = TT and can be lifted by breaking the inversion 
symmetry between the left and right domains. We fur- 
ther find that a local potential Uioc on B^^ or A^^ can 
raise and lower the energies of the green bands. Similarly, 
a line potential on B^^ or A§^ can change the energies 
of the magenta bands. In view of these results we pro- 
pose a criterion controlled by a hierarchy of energy scales 
to determine the number of fragile gapless modes in the 
atomically abrupt LSW shown in Fig. [5j 
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where lS.gap/'2. is half size of the field-induced gap which 
saturates if V exceeds a critical value [5]. Two gapless 



channels emerge if Eq. ( 12 ) is satisfied, but extra gapless 



channels can also appear from the flat bands if Eq. ( 13 ) 
is fulfilled. 

When the LSW is made smooth in the sense that it 
does not produce sufficiently strong intervalley coupling, 
the tunneling amplitudes near the domain wall in both 
layers are almost the same as the pristine ones. In such 
a case, the tunneling between the left and right domains 
in the broken (continuous) layer would strongly split the 
two green (magenta) bands at ka = tt. As a result, only 
one green and one magenta bands in Fig. [7] survive in the 
band gap, recovering our earlier continuum results. 

We conclude that the gapless interface modes at a LSW 
are topologically stable only if the potential difference be- 
tween layers is the dominant energy scale, so that valley 
is approximately a good quantum number. In the general 
case the number of domain wall modes can be any integer 
from to 4 depending on the criteria like that implied by 
Eq. (12) and (13). The valley-projected topological-state 



physics of BLG is illustrative of similar physics which 
occurs in all multi-layer graphene systems [3 [12] and is 
sensitive to stacking order, and to perpendicular electric 
fields. 

Note added. — After the finalization of this work, a 
complementary preprint [18j , which covers closely related 
material, has appeared. 

This work is supported by DARPA under grant 
SPAWAR N66001-11-1-4110, by the Department of En- 
ergy, Office of Basic Energy Sciences under contract 
DE-FG02-ER45118, and by the Welch Foundation grant 
TBF1473. 



6 



* |mele@physics.upenn.edu| 
[1] H. Min, A. H. MacDonald, Phys. Rev. B 77, 155416 
(2008). 

[2] E. McCann and V. I. Fal'ko, Phys. Rev. Lett 96, 086805 
(2006). 

[3] T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Roten- 

berg, Science 313, 951 (2006). 
[4] E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 

Peres, J. M. B. Lopes dos Santos, J. Nilsson, F. Guinea, 

A. K. Geim and A. H. Castro Neto, Phys. Rev. Lett. 99, 

216802 (2007). 

[5] F. Zhang, B. Sahu, H. Min and A. H. MacDonald, Phys. 

Rev. B 82, 035409 (2010). 
[6] B. E. Feldman, J. Martin and A. Yacoby, Nature Physics 

5, 889 (2009). 

[7] J. Martin, B. E. Feldman, R. T. Weitz, M. T. Allen and 
A. Yacoby, Phys. Rev. Lett. 105, 256806 (2010). 

[8] E. A. Henriksen and J. P. Eisenstein, Phys. Rev. B 82, 
041412(R), (2010). 



[9] J. Velasco, L.Jing, W. Bao, Y. Lee, P. Kratz, V. Aji, M. 
Bockrath, C. N. Lau, C. Varma, R. Stillwell, D. Smirnov, 
F. Zhang, J. Jung, and A. H. MacDonald, Nature Nano. 
7, 156 (2012). 

[10] W. Bao, J. Velasco, F. Zhang, L. Jing, B. Standley, D. 
Smirnov, M. Bockrath, A. H. MacDonald, and C. N. Lau, 
Proc. Nat. Acad. Sci. 109, 10802 (2012). 

[11] I. Martin, Y. M. Blanter, and A. F. Morpurgo, Phys. 
Rev. Lett. 100, 036804 (2008). 

[12] F. Zhang, J. Jung, G. A. Fiete, Q. Niu, and A.H. Mac- 
Donald, Phys. Rev. Lett. 106, 156801 (2011). 

[13] D. Xiao, M. C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 
1959 (2010). 

[14] J. Li, A. F. Morpurgo, M. Biittiker, and I. Martin, Phys. 

Rev. B 82, 245404 (2010). 
[15] J. Li, I. Martin, M. Buttiker, and A. F. Morpurgo, Nat. 

Phys. 7, 38 (2011). 
[16] Z. Qiao, J. Jung, Q. Niu, and A. H. MacDonald, Nano 

Lett. 11, 3453 (2011). 
[17] J. Jung, F. Zhang, Z. Qiao, and A. H. MacDonald, Phys. 

Rev. B 84, 075418 (2011). 
[18] A. Vaezi, Y. Liang, D. H. Ngai, L. Yang, and Eun-Ah 

Kim, |arXiv:1301. 1690 (2013). 



